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Abstract 

We calculate the cross-section for the excitation and subsequent decay of 
triple giant resonances (TGDR) in several nuclei excited with heavy ions. 
The recently developed coherent plus incoherent theory for the excitation 
in conjunction with the hybrid decay model of Dias-Hussein-Adhikari are 
used for the purpose. It is emphasized that the direct decay of the TGR 
is expected to deviate appreciably from the harmonic limit especially at 
low bombarding energies, owing to the incoherent contribution. 
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The study of the double giant dipole resonance in nuclei has received a considerable 
amount of attention over the last 15 years [1]. Both the pion double charge exchange 
and relativistic heavy ion Coulomb excitation reactions have been used to probe this 
large amplitude collective motion in may fermion systems. The quest for the similar 
double plasmon resonance in metallic clusters is underway [2] . Plans are also in progress 
to search for the triple giant dipole resonance (TGDR) in nuclei [3]. It is clearly of 
importance to supply theoretical estimates of the cross-section as well as the different 
decay branching ratios of these exotic collective modes. This is the purpose of the 
present paper. We use the recently developed coherent plus incoherent excitation 
theory of Ref. [4] in conjunction with the hybrid decay model of Dias-Hussein-Adhikari 
(DHA) of Ref. [5]. 

The existing models for the calculation of the excitation cross-section of DGDR 
can be grouped into three categories: a microscopic structure model in conjunction 
with second order Coulomb excitation perturbation theory [6], a macroscopic, oscilla- 
tor (harmonic or anharmonic) model in conjunction with coupled channels Coulomb 
excitation theory [7] and finally the recently developed average plus fluctuation model 
[4,8]. In this latter model the average cross-section is calculated according to the the- 
ory developed in [9], where the simple, double, etc. giant resonances are considered as 
doorway states belonging to the spectrum of a damped harmonic oscillator. 

A fully microscopic structure calculation of the excitation cross-section of the 
TGDR resonance is prohibitively difficult. A huge number of three particle-three hole 
configurations have to be dealt with in a coupled channel context. A detailed account 
of the spreading of the TGDR would require the inclusion in the calculation of at least 
the four particle - four hole subspace. Not having available such a detailed descrip- 
tion we opt for using our coherent + fluctuation model [8]. The TGDR excitation 
cross-section is found to have the form 

*V=*V+*f{2) + *%{l) (1) 
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where a c is the average cross-section for the coherent excitation of the three phonons 
which proceeds through the one and two phonon states in a typical three-step descrip- 
tion. This cross section can be guaranteed for the coupled channel Coulomb excitation 
model which contains explicit reference to flux loss from there states owing to their 
spreading into more complex configurations. The cross-section cr^ (2) corresponds to 
fluctuation contribution arising from the decay of an intermediate collective GR phonon 
into the complicated background followed by the excitation of another collective phonon 
on the background states (the Brink- Axel phonon). The final states involved in this 
cross-section contains two collective phonons in contrast to a c . Finally (1) con- 
tains contributions that lead to only one phonon in the final state. In terms of the 
time sequence of events, af \ being a three-step process, is the fastest, followed by 
the four-step process accounted for by crj^(2) (this is a four-step process since besides 
the three excitation steps one has one internal mixing step) and finally the five step 
process contained in cr^(l). The DGR cross section may be similarly decomposed as 



where Y\ is the spreading width of the single phonon GR and r c is the average collision 
time given by — , with b being the grazing impact parameter, v the asymptotic relative 



The above fluctuation contribution become insignificant at very high energies where 
the systems proceeds very quickly though the sequence ground state — > GDR 
DGDR — * TGDR. When the bombarding energy is lowered, the fluctuation effects 
may become appreciable [8]. The estimates above were found to be quite reasonable 





(2) 




(3) 
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when compared to the more elaborate model of Ref. [4] at relatively high bombarding 
energies. At lower energies, clearly, for the evaluation of the different contributions to 
the total excitation cross-section, Eq. (1), one should rely on the latter, more precise 
model. In order to discuss the decay of the final states into the open channels, we need 
to know the values of the different contributions to and a^ 3 \ 

We have calculated the excitation cross-sections, , and , for various nu- 
clei incident on 208 Pb at several bombarding energies, using a three-dimensional (3D) 
generalization of the model of Ref. [4]. The 3D time evolution equation used to de- 
scribe the excitation and decay of the GDR phonons possesses the same form as the 
one-dimensional equation of Ref. [4]. However, the collective and statistical excited 
states of the 3D model take into account all possible combinations of the (two) trans- 
verse and (one) longitudinal degrees of freedom, which yield 3 coherent one-phonon 
states, 6 coherent two-phonon states and 10 coherent three-phonon states as well as a 
multitude of states containing a mixture of coherent and statistical excitations. De- 
cays of the three types of phonons to the statistical background are assumed to occur 
independently but to each obey Bose-Einstein statistics. 

The Coulomb interaction matrix elements used to describe the transverse modes of 
the GDR excitation in the 3D model are the physically appropriate ones, as given in 
Ref. [9]. The longitudinal Coulomb interaction matrix element, however, is modified 
from the form given there. It is reduced to a term proportional to the longitudinal 
component of the eletric field, in analogy to the transerse terms, but which differs 
from the expression given in Ref. [9] by a total time derivative. The latter term can 
be extracted from the equations and discarded when only the coherent excitation is 
included. This is no longer the case when decay to the statistical states is taken into 
account. Nevertheless, we have neglected its contribution here. 

As in Ref. [4], the coupled equations of motion are solved as a function of impact 
parameter to yield asymptotic occupation probabilities. Effective asymptotic occupa- 
tion probabilities are defined, for states that decay, as the sum over the probability 
that decays out of each state during the time evolution. Cross sections are obtained 
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by integrating each probability x differential area over impact parameter and summing 
over polarizations. 

The various contributions to the cross sections are easily extracted from the theo- 
retical calculations. In Table I, we present the coherent and fluctuation contributions 
to the DGDR cross section, and cr^(l) for various nuclei incident on 208 Pb at 
several energies. In Table II, we present the contributions to the TGDR cross section, 
<jf \ 0-j^(2), and cr^(l). We observe that the cross sections increase dramatically with 
the charge of the projectile. As is well known, the coherent two-phonon cross sections 
scales approximately as the charge Z squared, while the three-phonon one scales as Z z . 
We also observe that the coherent contribution to the cross sections only dominates 
at relatively high incident energies. At the lowest energy shown, 10 MeV/nucleon, the 
cross sections involving statistical phonons are substantially larger than the coherent 
ones. At E/A — 100 MeV, it is clear from the tables that the fluctuation contribution 
to the DGDR cross section is about as large as the coherent one, while the fluctuation 
contribution to the TGDR is about three times larger than the coherent contribution. 

We turn now to the decay of the DGR and TGR. We first remind the reader of the 
hybrid direct +fluctuation decay model of DHA [5]. According to this model, which 
has been extensively used in the analysis of the decay data [10,11], the GR decays to 
a find channel / in the following manner: 



(GR) (CN) . (GR) 

— TTTrr; + AH" 



= a' 1 ' (P< + Pj) , (4) 

where cr ' is the one phonon excitation cross section discussed before, while fj,\ — ^ and 
the r's are the appropriate transmission coefficients. We have written the probability 
of populating the final channel / through direct decay of the GR as 

T (GR) 



5 



and the probability of of populating the channel / through the statistical states as 



Note that the statistical decay component contains explicit reference to the GR direct 



transmission, ^/i 1 rj G ' i? ' ) ^ 



Before entering into the details of the decay of the multiple giant resonances, let us 
first analyze the decomposition into direct decay and decay into the statistical states. 
For this purpose, we use the branching ratios 

P T = P f = 1 - Mi and pl = Y1 P f = Mi • ( 7 ) 

/ / 
The decomposition of the single GR decay is a direct result of Eq. (f|) , 

a W _ (pT + pi) a (D . ( 8 ) 

To decompose the decay of the multiple giant resonances into direct and statistical 
parts, we assume that each of the collective phonons decays independently. The decay 
of the coherent contributions to the DGR and TGR can then be decomposed as 

a f) ( P T + piy a {2) (9) 

= (P T2 + 2P T P i + P i2 )af , 



and 



cr. 



(3)^(pT + pl)3 (J (3) (1Q) 



(3) 



= (P T3 + 3P n P l + 3P T P i2 + P 13 ) o\ 

That is, the coherent contribution to the DGR can decay through direct decay of each 
of the collective phonons, through a direct decay of one of the colective phonons and 
decay into the statistical states of the other, or through decay into the statistical states 
of both of the phonons. Decay of the coherent contribution to the TGR takes into 
account the different possible direct or statistical decays of the three initial phonons. 
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We can analyze the decomposition of the decay of the fluctuating contributions to 
the DGR and TGR cross sections in a similar manner. We need only take into account 
the number of collective phonons in each of the contributions. Thus, the decay of the 
fluctuating component of the DGR can be decomposed as 

^W-^ + ^V/fW- (ii) 

while the decay of the fluctuating components of the TGR can be decomposed as 

af (2) - (P T2 + 2P T P^ + P 12 ) af}{2) , (12) 

and 

af(l)^(p1 + pl)vf(l). (13) 

We can now combine the various terms into decompositions of the decay of the 
total DGR and TGR cross sections, 

+ PT(2P'<7?> + <$(!)) 
+ pl(pivW+*%\l)), 

and 

a m^ P ^ a (3) (15) 

+ pT2 (3 p| (T (3) +(7 (3) (2)) 

+ P T (3 p| 2cr (3) +2 p V (3) (2)+(7 (3) (1)) 

+ pi(pi2 (7 (3) + pi (7 (3) (2)+(7 (3) (l))) 

where we have collected terms according to the number of direct decays involved. 

In medium to heavy nuclei, one expects the spreading to dominate over the escape 
from the GR, which implies P^ = /ii ~ 1. In this case, we see that the completely 
statistical DGR and TGR decay cross sections, given by the last term in the preceding 
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two equations, will be approximately proportional to their total excitation cross sec- 
tions. They will not distinguish between the coherent and fluctuating components of 
the cross sections. Further analysis of the decay of the statistical component to equilib- 
rium can be quite complicated. Usually, however, particle emission from the statistical 
component can be well described using the equilibrium Hauser-Feshbach formalism. 

To best view the distinction between the coherent and fluctuation contributions to 
the GDR cross sections, we look for effects in the direct decay from the giant resonance. 
We have already analyzed the effects of the coherent and fluctuating components on 



exclusive decays of the DGDR. | 12| . Here, we wish to analyze their effects on inclusive 
decay cross sections of the DGR and TGR. 

To separate direct contributions to an inclusive cross section from the statistical 
contributions, it is necessary to concentrate one's attention on the high-energy end of 
the emission spectrum. There, the statistical weight of the high energy residuals in 
the statistical cross section strongly suppresses emission, leaving the direct emission 
to dominate the spectrum. To obtain the inclusive emission cross section from our 
decomposition of the decay of the excitation cross sections through direct and statistical 
modes, Eqs. [14] and [TR we take into account that each of the direct decay factors P^ 
can contribute independently to direct emission in channel / with a relative probability 
Pj I ' P\ Thus, for example, each of the factors of P^ in the P^ 2 (first) term of Eq. [14] 
contributes to the direct emission in channel / with probability Pj /P\ resulting in a 
total contribution of 2Pj / P^ of the term to the emission cross section. In general, we 
include a factor of Pj / P^ for each of the factors of P^ appearing in the equations. We 
neglect the terms containing only statistical decays, since we have Pj/P^ w at the 
high end of the emission spectrum. 

We point out that the relative decay probabilities Pj / P 1 ' are not quit the same for 
each of the decays, since conservation laws constrain the energies and angular momenta 
of the emitted particles to values consistent with the residual values in the nucleus. The 
residual nuclear values will not be identical for a cold nucleus (first emission) and a hot 
one (second and later emissions). Nevertheless, the general structure of the phonon, 
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which we assume to be essentially a one particle-one hole state, assure that the value 
of Pj/P^ will be about the same for each of the emissions. 

Using the rules above for estimating the contribution of the P T factors to emission, 
we find the direct inclusive emission spectra to be given by 

af = 2^P^4 2) (16) 

T 



and 



which reduce to 



and 



af = Z?LpVaV> (17) 
+ ^pt(3P^Vf +2P^>(2) + <$(!)), 



a? = P} (2a™ +*$(!)), (18) 



af = P}^ + 2af{2) + af{l)). (19) 

Here we see clearly the importance of the (more) coherent contributions to the direct 
emission spectra. Each of the excitation cross sections contributes according to the 
number of collective phonons it possesses. 

We can similarly analyze the contribution to the inclusive cross sections of each of 
the components of the of the excitation cross sections. In particular, we determine the 
contribution of the coherent DGR and TGR to the emission spectra to be 

<r% = 2P}*V and af f = ZP}af\ (20) 

We could call these values the harmonic limit of the cross section. Comparing these to 
the values for the total direct emission, we find 

J— = 1 + - fl w (21) 

(2 ) i + 2 ( 2 ) > y Zi ) 



c 



and 




3 ai 3) 3 



2 (2) , la 




(22) 



Thus, a considerably larger direct decay may occur if the fluctuation contributions are 
important, which may occur at lower bombarding energies. Of course, one could ob- 
tain deviation of the direct decay from the harmonic limit (two or three independently 
decaying phonons), if anharmonic effects were allowed, This, however, will imply devi- 
ation of the spectrum of the oscillator from the harmonic sequence, which seems to be 
borne out neither by experiment [1] nor by calculation [7]. 

In conclusion, we have, in this paper, calculated the excitation cross-section and 
studied the the decay properties of the double and triple giant dipole resonances of 
various nuclei as excited in Coulomb collisions with 208 Pb. It was found that the 
degree of deviation of the direct decay from the limit of two or three independently 
decaying collective phonons depends significantly on the bombarding energy, E, and 
can be appreciable at low values of E. 
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TABLES 

TABLE I. Contributions of the coherent and fluctuation components to the DGDR 
excitation cross section (in mb) of various projectiles incident on a lead target at several 
different values of the incident energy. 



Projectile 


<£> 


10 MeV 


100 MeV 
4 2) <$(!) 


1 GeV 

aP 




40 Ca 


1.01 


4.52 


2.17 


2.19 


7.20 


0.72 


120 Sn 


12.37 


49.54 


26.48 


22.94 


72.61 


6.65 


132 Xe 


15.03 


59.62 


32.19 


27.57 


88.50 


8.00 


165 Ro 


23.32 


90.61 


51.13 


42.60 


138.59 


12.34 


208p b 


40.45 


148.66 


96.95 


72.87 


234.84 


19.83 


238 |j 


47.05 


174.78 


109.15 


84.86 


276.53 


24.04 



TABLE II. Contributions of the coherent and fluctuation components to the TGDR 
excitation cross section (in mb) of various projectiles incident on a lead target at several 
different values of the incident energy. 



Projectile 


aP 


10 MeV 

4? ( 2 ) 




aP 


100 MeV 

4? ( 2 ) 




aP 


1 GeV 

4? (2) 


4?« 


40 Ca 


0.01 


0.08 


0.11 


0.02 


0.06 


0.02 


0.11 


0.02 


0.00 


120 Sn 


0.17 


2.46 


3.17 


0.84 


1.92 


0.64 


3.03 


0.47 


0.04 


132 Xe 


0.23 


3.20 


4.10 


1.10 


2.50 


0.83 


4.07 


0.62 


0.05 


165 Ro 


0.42 


5.81 


7.37 


2.08 


4.70 


1.54 


7.76 


1.17 


0.09 


208p b 


0.95 


12.24 


14.83 


5.28 


10.78 


3.36 


16.68 


2.40 


0.18 


238 |j 


1.18 


15.09 


18.60 


6.13 


13.14 


4.22 


21.01 


3.14 


0.24 
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